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Abstract 

Transverse ionization cooling of muons is modeled as a Brownian motion 

of the muon beam as it traverses a Li or Be rod. A Langevin like equation 

is written for the free particle case (no external transverse magnetic field) 

and for the case of a harmonically bound beam in the presence of a focusing 

magnetic field. We demonstrate that the well known muon cooling equations 

for short-absorbers can be extrapolated to the useful case of a long absorber 

rod with a focusing magnetic field present. 
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I. INTRODUCTION 



The possibility of a collider to explore the Higgs energy range and supersymmetry 

has begun to be vigorously examined. One of the crucial issues to achieve the required 
luminosity (C ~ 10 34 cm _2 s _1 ) is the need to compress the phase space by means of muon 
cooling. A technique that has been shown to be very promising is ionization cooling. The 
introduction of the concept and the physics was first discussed by Skrinsky for a clear 
and comprehensive treatment we refer the reader to Neuffer's articles 0. 

The original derivation in ref. of the transverse cooling differential equation assumed 
a cooling system consisting of small alternating absorber and reaccelerator sections. Subse- 
quently, Palmer || and Fernow |4j have argued that the equation is valid for a single long 
absorber. Their argument is quite straightforward; from the definition of emittance it is 
possible to show that 



def .dE M 1 ( 7 /3) 



dz 1 dz 'P 2 E, 2el{z 



9 d <9 2 > „ 9 d <r 2 > 
< r 2 > h < 9 > 



dz dz 
e]_ (z) dz 

where we have used ^p- = -p^s^j^-, e ± is the muons' normalized transverse emittance, 
is the muon total energy, /3j_ is the beta function and < r 2 >, < 9 2 > are the square of the 
rms position and divergence of the beam due to multiple scattering. The first term in Eq.|l| 
reflects the energy loss (cooling) and the last three terms are produced by multiple scattering 
( heating) . 

Assuming a long cooling rod (Li or Be) the gaussian approximation is quite adequate, 
then 

< y 2 > = \ey 

<e 2 > = e 2 c z (2) 

where the projected scattering angle o = 9 c y/z = 13,6 I Me ^ \ff- and Lr is the radiation 



(3cp y Lj 

length. In this expression we have neglected logarithmic correction terms ||. There is also, 
the more accurate formula for 9 D @ 



with the characteristic angle \c = V 0A57 ^ e Y}l^J MI±lL the phenomenological parameters 
F = 0.98 and v = 2(1*- f) are cnosen to fit the experimental data (Cl represent the mean 
number of scatters in the medium). 

If we neglect from Eq.p] the third and fourth terms, then we can write the minimum 
achievable emittance as [0], 

N _ (13.6[MeV}) 2 P±(0) ,,v 
±lmm ~ 2{3mcY-t\ L * 

Palmer || has pointed out that the third and fourth terms in Eq.|l| need not be included 
because the multiple scattering medium (Li or Be rod) is immersed in a uniform transverse 
magnetic field, which prevents the beam from spreading laterally. This raises the question, 
how is the functional form of the particle distribution changed in position and angle due to 
the external magnetic field? We will examine this question in the next sections. 

We should also point out that the treatment considered here is also of interest for a 
number of other problems in accelerator physics. These include scattering of particles by 
residual gas in a synchrotron JF], scattering in the Inverse Cerenkov accelerators ||, plasma 



beat- wave accelerators and plasma lenses for future linear colliders 10] and the dynamics 



of space-charge dominated beams [11 



II. PARTICLE DISTRIBUTION WITHOUT A MAGNETIC FIELD 



This problem has been analyzed in detail by Rossi [[R] following the Fokker-Planck 
equation approach. Let W(y, 9, z; y a , 9 Q ) dy d6 represent the number of particles in the phase 
space element (y, y + dy; 6,9 + dO) after traversing a medium of thickness z and initial 
coordinates y(0) = yo, 9(0) = 9 Q . It satisfies the equation 

— - ( *- 0o) ^ + ih^ (5) 

with boundary conditions W(y, 9, z\ y a , # )| 2= o = S(y — y )5{6 — 9 Q ) and solution 



\/3 

W(y, 9, z : y a , 6 ) = — exp 



2 ({e-Q f 3(y-Y )(9-e ) + 3(y-Y ) r 



(6) 



9 2 \ z z 2 z 3 

which can be verified by direct substitution with O = 9 Q and Y Q = y Q + 9 D z. This result 
allows us to compute the emittance of the beam after traversing the cooling rod. After 
tedious gaussian integrations we obtain 

<y> =Y Q 

<9> =9 Q 

9 2 z 3 

<y 2 > =Yo+-f- 

<9 2 > = 2 + 9\z 

9 2 z 2 

<y9> = B Y + -JL- (7) 

Averaging over the initial coordinates and assuming gaussian distributions with standard 
deviations a yo and o~g , we obtain 

<y> = < # >= 



9 2 z 3 

+ 2z < y o o > + 

2 ^ _2 , n2. 



< y 2 > = alo + °l z 2 + 2^ < y o o > 4 1/1 



« 9 2 » = al + 9 2 c z 



9 2 z 2 

< ^ > = a e 2 c z+ < y o o > (8) 



The total emittance in the absence of a focusing field is 



9 4 z 4 z 3 

e ± (z) = t /ei(0) + + oi^- + £ c 2 < > + a^z (9) 



The terms proportional to 9 C are the contributions due to multiple scattering. 



III. PARTICLE DISTRIBUTION WITH A MAGNETIC FIELD 



Consider now the problem of a particle traversing a rod of material that has an axial 
current flowing through it; such a particle satisfies the equation of motion j% + K(z)y = 
where K(z) = eB a = u 2 , B is the azimuthal magnetic field and a is the radius of the rod. 



mcyf3a 
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K(z) is a function of z because of the energy loss, but for the simplicity of the arguments 
that follow, we neglect the energy change as the beam traverses the rod. 

A more complete treatment must take into account random accelerations of the particles 
due to scattering (i.e. stochastic changes in angle = 9). A correct equation of motion is 

where we denote with A(z) the random acceleration due to Coulomb scattering which excites 
betatron oscillations in the beam. This equation is formally a Langevin equation of a particle 
in an external field K(z)y (harmonic oscillator) where the frequency is a function of the time 
variable z. The main assumptions regarding the stochastic variable A(z), more precisely 
jz+ dz dz'A(z'), is that it is independent of y, that it varies extremely rapidly compared to 
the variations of the coordinates y and 9, and that it is Gaussian-distributed with a variance 
01 

Therefore, casting the muon cooling problem in stochastic terms, we first determine the 
particle distribution W(y,9, z;y 9 ), and then as before, we calculate the emittance from 
that function. 

The method for determining the distribution function uses standard techniques for solv- 
ing ordinary differential equations and is described in detail by Chandrasekhar fl3| . The 
result for the distribution function is, 



W(y, 9, z,u;y ,u ) 



2nJ{FG-H 2 } 



cxp 



]{FC !_ H2} (G(y - Y a f - 2H(y - Y ){9 - O ) + F{9 - Q ) 



2\ 

where the parameters F,G H are functions of the external focusing field; 

z , sin2u;z N 



'ir 



e. 
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2a; 2 v 2uz ' 
2 z sin2u;z 

G -° c 2 {l + ^z- ) 
H = 9 2 c ^(l-cos2uz) 



e 



c J2 



FG — H = ^-z 

U0 A 



1- 2 

UJZ 



(12) 



It can be shown that Eq.|TT] reproduces Eq.|6] in the limit u — > 0, and the probability den- 
sity W(y,9, z,uj;y ,9 ) satisfies a parabolic partial differential equation, the Fokker-Planck 
equation 



-97 = -^- 0o) ^ + u (y ~ Yo) ^o- + (13) 



As in the previous section we are interested in calculating the second moments of the dis- 
tribution and from those the emittance. We find that 

sin ujz 

< y > — y cos ujz + 9 Q z 

UJZ 

< 9 > = 9 a cos uz — y Q uj sin uz 

sincuz.n 9 2 , sin2u;2:, 

Y + z—^(l 

uz ' 2u 2K 2ujz 



<y 2 > =(yo^z + 9 z—y + z^{i-^ r ) 

< 9 > = [9 cos ojz — y u sin ojz) + z— (1 + 



2 v 2cj2 



,2 



y"u; „ 9 sin ojz # 2 ,2 2 .sinu;;?. 9 , 

< ^ > = y 9 cos2ujz- Z°—sm2uz + (Pzcoawz h -£—{ ) (14) 

2 uz 2 ujz 

If we now assume an uncorrelated ensemble of incident particles with independent gaussian 

distributions of initial conditions y Q and 9 Q and average over both variables, we obtain 

< y > = < 9 >= 

9 9 . , 9 , ,.sinwz,, 91 ,„ sin2a;z. 

« r » = «4(o» w *)' + ) 2 + - - as-) 

< # 2 > = a 2 8o (cosuJz) 2 + a 2 uj 2 (sin ujz) 2 + z-f (1 + 

J 2 2w2 



» sin 2^ 2 1 . n 9 2 r z 2 .sinujz , 2 

« " » = ^-sr " sin2 " 2 + > (15) 



An important observation is that for a very high magnetic field (uj » 0) the rms beam 
size remains approximately constant. If the beam is focused to a waist at the entrance to 



the rod, the emittance at any distance z inside the rod is (15 



el(z) = ei(0) + a 2 yo 9 2 c z + ^z 2 + cos ( 2 ^)] ( 16 ) 

The second term gives the dominant contribution of multiple scattering to the emittance. 
This also leads to Eq.f| for the minimum emittance, provided that a 2 Q 3> ^2L rod and 
°yo ^ which will usually be the case for strong focusing (cu 0). 
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IV. CONCLUSIONS 



Using an analogy with a random dynamical process modelled with a Langevin equation, 
we have incorporated the stochastic nature of both position and angle variables into the 
problem of a muon traversing a Li or Be rod immersed in a uniform azimuthal magnetic 
field. The pseudo-Brownian motion of the particles in the medium represents heating of the 
beam. The emittance increase due to multiple Coulomb scattering (the less likely single and 
plural scattering events are neglected) opposes the emittance decrease (cooling), introduced 
by the energy loss With sufficiently strong focusing present, ionization cooling can 
effectively take place over long lengths of absorbing material. 
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